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INTRODUCTION 



There is by now convincing evidence that the available 
high quality cosmological data (Type la supernovae, 
CMB, etc.) are well fitted by an emerging "standard 
model". In the context of GR this "standard model", 
assuming flatness, is described by the Friedman equa- 
tion 

H 2 (a) = = Hq [n m (a) + Q DB (a)] (1) 

where a(t) is the scale factor of the universe, f2 m (o) 
is the density parameter corresponding to the sum of 
baryonic and cold dark matter, with the latter needed 
to explain clustering, and an extra component Ode (a) 
with negative pressure called dark energy needed to ex- 
plain the observed accelerated cosmic expansion (Riess 
et al. 1998; Perlmutter et al. 1999; Efstathiou et al. 
2002; Tegmark et al. 2004; Spergel et al. 2007; Nesseris, 
& Perivolaropoulos 2005; Nesseris, & Perivolaropoulos 
2007a). During the last decade there have been many 
theoretical speculations regarding the nature of the ex- 
otic "dark energy". Various candidates have been pro- 
posed in the literature, among which a dynamical scalar 
field acting as vacuum energy (Ozer & Taha 1987; Cald- 
well, Dave & Steinhardt 1998; Peebles & Ratra 2003). 
Under this framework, high energy field theories gener- 
ically indicate that the equation of state of such a dark 
energy is a function of the cosmic time. To identify this 
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type of evolution of the equation of state, a detailed 
form of the observed H(z) is required which may be ob- 
tained by a combination of multiple dark energy probes. 
Such probes may be divided in two classes according to 
the methods used to obtain H(z). 



• Geometric methods probe the large scale geometry 
of space-time directly through the redshift dependence 
of cosmological distances [(1l(z) or cIa(z)]. They thus 
determine H(z) independent of the validity of Einstein 
equations. 

• Dynamical methods determine H(z) by measuring 
the evolution of energy density (background or pertur- 
bations) and using a gravity theory to relate them with 
geometry ie with H(z). These methods rely on knowl- 
edge of the dynamical equations that connect geometry 
with energy and may therefore be used in combination 
with geometric methods to test these dynamical equa- 
tions. 

A very accurate and deep geometrical probe of dark 
energy is the angular scale of the sound horizon at the 
last scattering surface as encoded in the location l± T 
of the first peak of the Cosmic Microwave Background 
(CMB) temperature perturbation spectrum. This probe 
is described by the so called CMB shift parameter 
(cf. Bond, Efstathiou & Tegmark 1997; Trotta 2004; 
Nesseris & Perivolaropoulos 2007) which is defined as 
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Figure 1. The locus of v)q, W\ for Case 1 (solid line), Case 2 (dotted line), Case 3 (circle, see text), Case 4 (m = 2, n = short 
dashed line) and Case 4 (m = 1, n = 1 long dashed line). 
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where if T is the temperature perturbation CMB spec- 
trum multipole of the first acoustic peak. In the defini- 
tion of R, lJ T corresponds to the model (with fixed fi 
fib and h) characterized by the shift parameter and Z x 
to a reference fiat SCDM model (f2 m = 1) with the same 
uj m — Qmh 2 and cj b — Q, b h 2 as the original model. Re- 
cently, Nesseris & Perivolaropoulos (2007) have found 
that models based on general relativity that have iden- 
tical shift parameter R and matter density f2 m also lead 
to almost identical ISW effect despite of their possible 
differences in the cosmic expansion histories. 

The aim of this work is to investigate our suspicion 
that the (geometrical) CMB shift parameter is somehow 
associated with the (dynamical) fluctuation growth rate 
in the context of general relativity. Note, that a possible 
violation of this connection may thus be viewed as a hint 
for modifications of general relativity. The structure of 
the paper is as follows. The basic theoretical elements 
are presented in section 2. The results are presented 
in section 3 by solving numerically the time evolution 
equation for the mass density contrast for various flat 
dark energy models that share the same value of shift 
parameter and value of Q m . In section 4 we draw our 
conclusions. Finally, in the appendix we have treated an- 
alytically, up to a certain point, the differential equation 
for the mass density contrast considering different dark 
energy models with a time varying equation of state. 



2 THEORETICAL ELEMENTS 

The location T of the first acoustic peak in the CMB 
temperature spectrum can be connected with the angu- 
lar diameter distance d,A and with the sound horizon r s 
both at the last scattering surface (z = zi 8 ) and then the 
shift parameter (see Nesseris & Perivolaropoulos 2007 
and references therein for details) can be brought to the 

form R' = 2 - . The expression usually used 
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for the shift parameter is 
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where E{z) = H(z)/H . 

In order to define E(z) we use the Chevalier- 
Polarski-Linder (CPL) [Chevallier & Polarski 2001; Lin- 
der 2003] parametrization for which 



w(a) = Wo + Wi(l — a) 



and 
E\a) 
where 
f(a) = exp 



Q m a- 3 + (1 - fim)a " 



/(a) 
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On the other hand, a dynamical probe of geometry 
is the measured linear growth factor of the matter den- 
sity perturbations 5(a). The evolution equation of the 
growth factor for models where the dark energy fluid 
has a vanishing anisotropic stress and the matter fluid 
is not coupled to other matter species is given by (Pee- 
bles 2003; Stabenau & Jain 2006; Uzan 2007): 
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Figure 2. The growth factor for the first three cases and all Monte-Carlo values of wq and w\ (fi m = 0.25). Case 1 corresponds 
to Fig. 2a and cases 2 and 3 to Fig.2b and Fig.2c respectively. Note that for this figure we do not normalize the growth factor 
with its value at a = 1 as we are interested in its range of values at z = 0. 



5 + 2H(t)S - AnGpmS = 
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where dots denote derivatives with respect to time. Use- 
ful expressions of the growth factor can be found for 
the ACDM cosmology in Peebles (1993) and for the 
quintessence scenario (to — const) in Silveira & Waga 
(1994), Wang & Steinhardt (1998), Basilakos (2003), 
Nesseris & Perivolaropoulos (2008) and for the scalar 
tensor models in Gannouji & Polarski (2008). As an 
example, in the case of ACDM cosmology the growth 
factor is of the form 
5fi n __ 

3 ) 3 
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which for a — 1 has some similarity with the form of 
the shift parameter (see eq.3). Finally, it is interesting 
to mention here that for dark energy models with a 
time varying equation of state, an analytical solution for 
5(a) has yet to be found because the basic differential 
equation (7) becomes more complicated than in models 
with constant w. However, in a recent paper (Linder & 
Cahn 2007) a growth index 7 was used to parameterize 
the linear growing mode including models with a time 
varying equation of state (see next section). 



3 THE CONNECTION BETWEEN THE 
CMB SHIFT PARAMETER AND THE 
GROWTH FACTOR 

In order to explore the above mentioned connection of 
the CMB shift parameter R of eq. (3) with the growth 
factor 5(a), we perform a Monte-Carlo analysis in the 
parameter space of wo — wi of the parametrization (4). 
In particular we fix fi m (for example 0.25) and compare 
the variation of the growth factor for wo — wi pairs cor- 
responding to fixed R with the corresponding variation 
when other combinations of w — wi are fixed. Specifi- 
cally, we consider and compare 4 cases: 



• Case 1: wo and wi are constrained by a fixed value 
of the CMB shift parameter, R(wo,wi) = 1.7 

• Case 2: wo and wi are constrained by a linear re- 
lation of the form wi = —2.9 — 3.1wo approximating 
the locus of the wo and wi that satisfy the relation 
R(wo,wi) = 1.7 (see Fig.l). This approximation is ac- 
curate to within about 5% and provides an estimate 
of the uncertainties introduced in the predicted value of 
the growth factor if the shift parameter is not accurately 
measured. 

• Case 3: wo and wi are constrained to be on a circle 
of radius 0.5 and center the ACDM point (-1,0) 

• Case 4: wo and wi are constrained by a fixed value 
of an integral ansatz of a form similar to the CMB shift 
parameter 

dz 



A(w , wi) r> 



(9) 



, (J (1 + z) n E m (z) 

for various values of the parameters m,n (see Table 1). 
Note that for (m, n) = (1, 0) we get the usual CMB shift 
parameter. 

In the first two cases wo is a random variable in the 
range [—1.5,-0.5] and for each value of wo we use the 
constraint equation to solve for wi. In case 3, wo and 
wi are on a circle of radius 0.5 and center the ACDM 
point (-1,0). In Fig.l we show the locus of the points 
wo and wi that satisfy the constraints of cases 1-4. No- 
tice that in Fig. 1, case 3 is shown as an ellipse and 
not a circle due to the fact that the aspect ratio is 
chosen so that the loci for the other cases are shown 
optimally. Also, in Fig. 2 we present the growth factor 
evolution which is derived by solving numerically eq. 
(7), for the first three cases and all Monte-Carlo values 
of wo and wi (100 pairs each time). In case 1 the values 
of 5(a = 1) are much more constrained than the other 
two cases as the range (dispersion) of the growth factor 
at a = 1, ie Max[S(a = 1)] — Min[5(a = 1)] for the 
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first case is 0.055 (corresponding to a variation of the 
mean value of 5(a = 1) of less than 8%) while in the 
second and third cases it is 3.1 (33% variation) and 5.9 
(49% variation) times that value. We have also investi- 
gated the sensitivity of our analysis to the matter den- 
sity parameter. In particular, we confirmed that in the 
range fi m G [0.2, 0.3] our results depend weakly on the 
value of fi m . In fact, the present time dispersion of the 
growth factor for fixed shift parameter varies from 3.5% 
for fi m = 0.2 to 9% for !T2 m = 0.3 (8% for fi m = 0.25). 
Thus our main result persists for all physical values of 
flm and it strongly indicates that the CMB shift param- 
eter is somehow associated with the growth factor in the 
context of general relativity. 

The linear relation of the form wi = —2.9 — 3.1u>o 
corresponding to Fig. 2b provides a rough approxima- 
tion (good to about 5%) of the locus of points that sat- 
isfy the relation R(wo,wi) = 1.7. This introduces sig- 
nificant additional dispersion to the present day growth 
factor (the dispersion goes to 33% in Fig. 2b from the 
8% obtained with the exact locus of fixed R in Fig. 2a). 
Once we improve the wi expression with an appropriate 
quadratic term in wo, the approximation improves from 
about 5% to about 0.3% and the growth factor disper- 
sion drops back to 8.5% (almost the same as with the 
exact locus of fixed 7?). This is an interesting result as 
it means that the shift parameter should be measured 
with la errors better than 1% for a determination of 
the growth factor to an accuracy better than 10%. This 
accuracy of measurement of the shift parameter how- 
ever is not far from present day measurements which 
have determined R to within 1.5% (R = 1.7 ±0.03 from 
Wang & Mukherjee 2007). 

In case 4 the analysis is similar to case 1, ie wo is a 
random variable but now wi is found from the general- 
ized constraint 



A(w ,wi) ■r 



A 



fixed m,n 



(10) 



where Af ixed m>n is the value of A(wo, wi) m ,n for a 
ACDM cosmology (wo = — l,wi = 0) for the respec- 
tive values of m, n. We computed the range (disper- 
sion) of the growth factor at a = 1, ie Max[S(a — 
1)] — Min[8(a = 1)] for the values wo and wi derived 
from the constraint eq. (10). In Table 1 we show the 
ratio of these values for various m, n to the value of the 
range of the growth factor for m = l,n = [when 
A(wo,wi) m ,n goes over to the CMB shift parameter 
R],ie 



(Max[5{a = 1)] - Min[8(a = 1)]) ( 



Table 1. The ratio of the values of the range of the growth 
factor at a = 1 for various m,n to that for m = l,n = 0. 
The case m = 1, n = corresponds to the CMB shift param- 
eter. Notice that for these values of m, n the ratio becomes 
minimal. 
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9.64 


9.74 


9.75 


9.16 


9.79 


1 


1.00 


5.92 


10.71 


11.72 


54.24 


2 


7.95 


10.84 


11.58 


56.81 


376.47 


3 


10.89 


11.89 


69.62 


428.33 


2969.39 


4 


11.61 


73.53 


247.44 


1422.99 


20633.90 



Finally, we checked that random values for both wo 
and wi give a much larger dispersion on the values of 
5(a = 1), as exactly they should. Notice that the CMB 
shift parameter R (m = \,n = 0) seems to constrain 
the growth factor by about an order of magnitude or 
more compared to other forms of the integral ansatz 
(different values of m and n). From a theoretical point 
of view a possible relation between the CMB shift pa- 
rameter and the growth factor can be used as a viable 
test for the general relativity. As expected (Bertschinger 
2006), changing the validity of Einstein's field equations 
(the so called theory of modified gravity), we change 
accordingly the growth factor [see Gannouji & Polarski 
(2008) for a detailed investigation of the growth factor 
in scalar-tensor theories]. In contrast the behavior of the 
CMB shift parameter remains unaltered, simply because 
the latter is a geometrical function (see e.g. Nesseris & 
Perivolaropoulos 2007). Thus, a mismatch between the 
measured value of the shift parameter and the measured 
value of the linear growth factor would be a hint towards 
modified gravity. 

Verifying this connection between the growth factor 
and the CMB shift parameter R analytically requires an 
exact or approximate solution to the differential equa- 
tion for the evolution of density perturbations when 
evolving dark energy is taken into account. In what fol- 
lows we discuss some recent attempts towards the con- 
struction of such approximate solutions. 

A well known approximate solution to eq.(7) is 
found by Linder & Cahn (2007), where a growth index 
7 was used to parameterize the linear growing mode for 
models with a time varying equation of state. Specifi- 
cally, the growth index 7 was defined through 

fiZ.(«) 



D(a) 



exp 



-du 



(12) 



where D(a) is the growth factor normalized to unity and 
it was found that 7 can be approximated by 
6-3(l + Woo ) 



m.n) 



(Max[S(a = 1)] - Min[S(a = 1)]) (1>I 



(11) 



7 11-6(1 + ^00) 
where ui-o = w(z » 1). 



(13) 
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Figure 3. Left panel: The growth factor as a function of redshift. The solid line represents our analytical approximation (see 
(14)), while the dashed line represents the parametrized growth factor (see (12) derived by Linder & Cahn (2007)). Right panel: 
The percent accuracy (1 — D appr / D num )% of the two approximations, cq. (12) (dashed line) and eq. (14) (solid line). Note, 
that we use f2 m = 0.25 and (u>o,u>i) = (—0.95,0.43). 



The previous approach was based on the approx- 
imation that the universe is not too far from being 
matter-dominated. However, by utilizing just some basic 
elements from the differential equation theory we have 
solved analytically up to a certain point eq.(7) (see the 
Appendix for details), assuming that the equation of 
state parameter is a function of time. In this approach 
the growing D + (a) and decaying D_ (a) modes are given 
by 

D±(o) * a-^E-^(a) exp (r^f J" | 3 (.)| 1/2 d^ (14) 

where g(a) is denned in the Appendix. 

Utilizing the best-fit cosmological parameters 57 m = 
0.25 and (wo,u>i) = (—0.95,0.43) obtained from the 
Gold06 SNIa dataset and the CMB shift parameter, in 
the right panel of Fig.3, we present the approximated 
growth factor (left panel), as a function of redshift by 
utilizing eq. (14) (solid line) and eq. (12) respectively 
(dashed line) . It is obvious that our analytical approxi- 
mation is indeed close to that found by Linder & Cahn 
(2007). In the right panel of Fig. 3 we show the percent 
accuracy (1 — D appr / D num )% of the two approximations 
Dappr, compared to the numerical solution. Eq. (12) de- 
viates from the numerical solution by 0.17% while our 
D+(z) approximation deviates by 0.5 — 1.1%. 



4 CONCLUSIONS 

In this work, we found that flat models which have the 
same fi m and CMB shift parameter R but different H(a) 
and w(a) also have very similar growth of perturbations 
even though the dark energy density evolution is quite 
different. This was done by comparing various forms of 



constraints for u>o and wi , besides the CMB shift param- 
eter, by a Monte-Carlo simulation. In all cases consid- 
ered, models constrained by the CMB shift parameter 
had also a very similar growth factor. 



APPENDIX 

In this appendix we try to treat analytically, as much 
as possible, the problem of the growth factor evolution 
for dark energy models with a time varying equation of 
state. 

The time evolution equation for the mass density 
contrast, modeled as a pressureless fluid, is obtained 
from the Euler and matter stress energy conservation 
equations as: 



S + 2H(t)S - AnGpmS = 



(15) 



where dots denote derivatives with respect to time. This 
differential equation is valid for models where the dark 
energy fluid has a vanishing anisotropic stress and the 
matter fluid is not coupled to other matter species, how- 
ever see Uzan (2007) for a detailed discussion of the 
modifications that appear on the right-hand side of the 
above equation when such terms are present. Changing 
variables from t to o the above equation becomes (see 
also Linder 2003) 



8 + A(a)8 + B{a)S = 
where A(a) = ^ 



(16) 



3 



i-n m +n m /(a) 



(l-n m )iu(q) 
[l-f2 m +n m /(a)J 



and B(a) 



Performing now the following transformation 



5(a) = y(a) exp 



A(u)du 



(17) 
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the linear mass density fluctuations is written: 

S(a) =y(a)a-' A/2 E- 1/2 (a) (18) 

and the unknown function y satisfies the following dif- 
ferential equation: 

y"-g(a)y = (19) 
with a relevant factor of 

g(a) = l -A{a) + \A 2 {a)- B(a) . (20) 

It becomes evident, that a major part of the pure solu- 
tion is described by the expression a~ 3 / 2 E^ 1 ^ 2 (a). 

Of course, in order to solve fully the problem we 
have to derive the functional form of y. In particular, 
we write eq. (19) as follows 



(21) 



\g(a)y 2 \ = \c + 



)du+y 2 \ (c 2 =2ci) .(22) 



From a mathematical point of view we can select the 
integration constant c to be large enough such as 

n\g(a)\y 2 <n 2 c 2 + [ n 2 \g(u)\y 2 (u) l -f^du+y 2 (23) 

where n 2 is the normalization constant of the prob- 
lem. Now, we can use Gronwall's theorem (see Gronwall 
1919), which is a well known theorem from the differen- 
tial equation theory. 

Theorem: Lets assume that fi : [a, f3] — > [0, oo) and 
y : [a, 0\ — » [0, oo) continuous functions and A € 71. If 

y<\+\ ( f i{x)y(x)dx\ Vte [a, /3] 

J t„ 

then 



(24) 



y < A exp ^| fi(x)dx\j \/t G [a, f3] (25) 
and using it on (23) we get the following useful formula: 
n 2 \g(a)\y 2 < n 2 C 2 exp (J° J§£|du) +y 2 . (26) 



Doing so it turns out that a possible approximation 
could be: 



■n\g{a)\y 2 ~ nc\g(a)\ +y 2 
from which we get that: 

y(a) ~ c cosh J |(/(it)| 1 ^ 2 du 

Taking into account eq. (18), we can obtain the following 
approximation: 



(27) 



(28) 



8{a) ~ c a- i,2 E- L ' 2 {a) cosh ( ±n j \g(u)\ 1/2 du ) .(29) 



In order to normalize our analytical expression we 
use as a limiting case the Einstein de-Sitter model 
(fi m = 1, g(a) = 21/16a 2 ) in which the behavior of the 
corresponding growing mode is well known D+(a) = a. 
Indeed, doing so we get n — \/2T/3 and thus, the fol- 
lowing normalized growing D + and decaying D_ modes 
respectively become: 

At(a) ~ a-^ 2 E-" 2 {a) exp (r^ J° Is^^^d^ (30) 
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